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The transverse-momentum distribution of the Higgs boson at the LHC 


The search for the Higgs boson is one of the highest priorities of the LHC physics program [jl|]. 
In the last years a significant effort has been devoted to refining the theoretical predictions for the 
various Higgs production channels and the corresponding backgrounds, which are now known to 
next-to-leading order accuracy (NLO) in most of the cases [2]. In the case of gluon-gluon fusion, 
which is the main Standard Model Higgs production channel, even next-to-next-to leading order 
(NNLO) QCD corrections have been computed, although in the large-M, approximation (M, being 
the mass of the top quark). The result has been obtained first for the total rate [[|], and more recently 
for fully exclusive distributions [3]. Among the possible observables, an important role is played 
by the transverse-momentum spectrum of the Higgs boson, whose knowledge may help to enhance 
the statistical significance of the signal over the background. 

When the transverse momentum qj of the Higgs boson is of the order of its mass Mu, the 
perturbative series is controlled by a small expansion parameter, as(Mjj), and the fixed-order pre¬ 
diction is reliable. The leading order (LO) calculation [t5| shows that the large-M, approximation 
works well as long as both Mu and qj are smaller than M,. In this framework, the NLO QCD 
corrections have been known for some time SUB- 

The small-gr region (qj <C Mu) is the most important, because it is here that the bulk of events 
is expected. In this region the coefficients of the perturbative series in a s (Mu) are enhanced by 
powers of large logarithmic terms. In m (Mj 1 /qf). To obtain reliable perturbative predictions, these 
terms have to be systematically resummed to all orders in as 101 • In the case of the Higgs boson, the 
resummation has been explicitly worked out at leading logarithmic (LL), next-to-leading logarith¬ 
mic (NLL) [10], [11] and next-to-next-to-leading logarithmic (NNLL) [12] level. The fixed-order 
and resummed approaches then have to be consistently matched at intermediate values of qj, so as 
to avoid double counting. 

In the following we present predictions for the Higgs boson q T distribution at the LHC within 
the formalism of Refs. [ 1_3]—[ L5]. In particular, we include the best perturbative information that is 
available at present: NNLL resummation at small qT and NLO calculation at large qj. 

An important feature of our formalism is that a unitarity constraint on the total cross section 
is automatically enforced, such that the integral of the spectrum reproduces the known fixed-order 
results. More details are given in Ref. [0. Other phenomenological results can be found in 
Ref. 


We now present quantitative results from Ref. [15] at NLL+LO and NNLL+NLO accuracy. At 
NLL+LO (NNLL+NLO) accuracy the NLL (NNLL) resummed result is matched to the LO (NLO) 
perturbative calculation valid at large qj. Our calculation is implemented in the numerical program 
HqT, which can be downloaded from ][F7 ]. The code is an improved version of the original program 
used in Ref. [LI], the main difference being in the matching procedure, which is now performed 
using the results of Ref. [$]. 

The numerical results in Figs. 1 and 2 are obtained by choosing Mu = 125 GeV and using the 
MRST2004 set of parton distributions [18]. At NLL+LO, NLO parton densities and 2-loop as are 
used, whereas at NNLL+NLO we use NNLO parton densities and 3-loop as. The NLL+LO results 
at the LHC are shown in Fig. [TJ In the left panel, the full NLL+LO result (solid line) is compared 
with the LO one (dashed line) at the default scales jlp = /Ir = Mu- We see that the LO calculation 
diverges to +°o as qj —> 0. The finite component, obtained through the matching procedure, is 
also shown (dotted line). The effect of the resummation, relevant below qj ~ 100 GeV, leads to a 
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Figure 1: LHC results at NLL+LO accuracy. 


physically well defined distribution at q T —> 0. In the right panel we show the NLL+LO band ob¬ 
tained by varying /.//.- and j.tp simultaneously and independently between 0.5 Mh and 2 Mh, imposing 
0.5 < f-F/fi-R < 2. The integral of the spectrum agrees with the total NLO cross section to better 
than 1%. The corresponding NNLL+NLO results are shown in Fig. |j. In the left panel, the full 



Figure 2: LHC results at NNLL+NLO accuracy. 

result (solid line) is compared with the NLO one (dashed line) at the default scales fLp = \Xr = Mh- 
The NLO result diverges to — °o as qp —> 0 and, at small values of qp, it has an unphysical peak that 
is produced by the numerical compensation of negative leading and positive subleading logarith¬ 
mic contributions. The finite component (dotted line) vanishes smoothly as qp —> 0, showing the 
quality of our matching procedure. The NNLL+NLO resummed result is slightly harder than the 
NLL+LO one, and its integral is in very good agreement with the NNLO total cross section. The 
right panel of Fig. |2| shows the scale dependence computed as in Fig. [TJ Comparing Figs. 1 and 2, 
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we see that the NNLL+NLO band is smaller than the NLL+LO one and overlaps with the latter at 
qj < 100 GeV. This suggests a good convergence of the resummed perturbative expansion. Other 
sources of perturbative uncertainty give smaller effects 101 - 

In summary, considering the above results, the perturbative uncertainty of the NNLL+NLO 
spectrum is of about 10 % at intermediate and small qj, where the bulk of the events is concentrated. 
At very small c /7 ( q T <10 GeV) non-perturbative effects should be taken into account, whereas at 
large c /7 the perturbative uncertainty increases. Our results for the qT spectrum are thus fully 
consistent with those on the total NNLO cross section [3j]. 
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